home *** CD-ROM | disk | FTP | other *** search
/ SGI Freeware 1999 August / SGI Freeware 1999 August.iso / dist / fw_plotutils.idb / usr / freeware / share / ode / coupled.ode.z / coupled.ode
Encoding:
Text File  |  1998-10-28  |  864 b   |  40 lines
  1. # This example simulates two similar pendulums, coupled with a horizontal
  2. # spring.  The motion is taken to be undamped.
  3. # 
  4. # The plot shows the displacement of one of the pendulums, as a function of
  5. # time.  It illustrates how energy (e.g. maximum displacement) oscillates
  6. # back and forth between the two.
  7.  
  8. # You may run this example by doing:
  9. #
  10. #    ode < coupled.ode | graph -T X -C
  11. # 
  12. # or alternatively, to get a real-time plot,
  13. #
  14. #    ode < coupled.ode | graph -T X -C -x 0 50 -y -0.6 0.6
  15.  
  16. # The equations are:
  17. # m*x1'' = - m*g*x1/l + k(x2 - x1)
  18. # m*x2'' = - m*g*x2/l - k(x2 - x1)
  19.  
  20. xone' = vxone
  21. xtwo' = vxtwo
  22.  
  23. vxone' = -g*xone/lone + k/mone * ( xtwo - xone)
  24. vxtwo' = -g*xtwo/ltwo - k/ltwo * ( xtwo - xone)
  25.  
  26. k = 1
  27. g = 9.8        # acceleration due to gravity
  28. lone = 5
  29. ltwo = 5
  30. mone = 1
  31. mtwo = 1
  32.  
  33. xone = 0.0
  34. xtwo = 0.3
  35. vxone = 0
  36. vxtwo = 0
  37.  
  38. print t, xone
  39. step 0,50
  40.